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ABSTRACT 

An on-going effort in the characterization of exoplanetary systems is the accurate determination of 
host star properties. This effort extends to the relatively bright host stars of planets discovered with 
the radial velocity method. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) 
is aiding in these efforts as part of its observational campaign for exoplanet host stars. One of the 
first known systems is that of 70 Virginis, which harbors a jovian planet in an eccentric orbit. Here 
we present a complete characterization of this system with a compilation of TERMS photometry, 
spectroscopy, and interferometry. We provide fundamental properties of the host star through direct 
interferometric measurements of the radius (1.5% uncertainty) and through spectroscopic analysis. 

We combined 59 new Keck HIRES radial velocity measurements with the 169 previously published 
from the ELODIE, Hamilton, and HIRES spectrographs, to calculate a refined orbital solution and 
construct a transit ephemeris for the planet. These newly determined system characteristics are 
used to describe the Habitable Zone of the system with a discussion of possible additional planets 
and related stability simulations. Finally, we present 19 years of precision robotic photometry that 
constrain stellar activity and rule out central planetary transits for a Jupiter-radius planet at the 5cr 
level, with reduced significance down to an impact parameter of 6 = 0.95. 

Subject headinqs: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (70 Vir) 


1. INTRODUCTION 

The exoplanet discoveries over the past couple of 
decades have revealed a particular need to understand 
the properties of the host stars. This is because the plan¬ 
etary parameters derived from the detection methods of 
radial velocities (RV) and transits rely heavily upon the 
mass and radius determinations of their parent stars. 
These are often determined using stellar models, but 
there are ongoing efforts to provide more direct measure- 
ments of the stellar properties through a steroseismology 
(iHuber et a .Il2014fl and interferom etry (jBovaiian et al.l 
I20I21120131 lyon Braun et 11111201411 . The importance of 
these measurements cannot be oyerstated since they not 
only affect the deriyed planetary properti es but also the 
quan t ification of the Habitable Zone (HZ) (|Kasting et al.l 
1 19931 iKopparapu et al.ll2013l 1201411 and subsequent cal- 
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culations of the fractio n of stars with Earth-size plan- 
ets in the HZ, or r?m jD r essing fc Charbonneaul 120131 
IKoppara^l2013t iPetigura et alji2013[) . 

The Transit Ephemeris Refinement and Monitoring 
Survey (TERMS) is aiding in stellar characterization for 
the brightest stars as part of it s program to improve ex- 
oplanetary orbita l parameters (|Kane fc von Braiinll2008l 
iKane et ^120091 1. Recent results include detailed spec¬ 
troscopic and phot ometric analyses of the planet-hosting 
stars HD 38529 (iHenrv et al.l [Ml and HD 192263 
(|Dragomir et al.l l2012f) and the identification of long¬ 
term activity cycles. TERMS observations have also 
led to the discove ry of new pla nets in the HD 3 7605 
(jWang et al.ll2012t l and HD 4203 (jKane et al.ll2014[ l sys¬ 
tems. These efforts are continuing with a focus on the 
brightest host stars, which tend to be those around which 
planets were discoyered using the RV method. 

One of the earliest exoplanet discoyeries was that of 
the planet orbiting the bright (V = 5) star 70 Virgi- 
nis (hereafter 70 Vir) . The planet was discovered by 
IMarcv fc Butled (1199611 and lies in an eccentric (e = 0.4) 
116 day orbit. iPerrvman et al.l l|1996ll subsequently used 
Hipparcos astrometry to constrain the inclination and 
thus determine that the companion is indeed sub-stellar 
in mass. Observations of 70 Vir have continued since 
discovery, with RV da ta from Observat oire de Haute- 
Provence published bv iNaef et al.l (12004( 1 and the com¬ 
plete L ick Observatory dataset compiled bv iButler et al.l 
(|2006f) . A detailed characterization of the nearest and 
brightest exoplanet host stars is important because these 
continue to be those which are most suitable for poten¬ 
tially studying exoplanetary atmospheres for transiting 
planets. 

Here we provide a detailed analysis of the 70 Vir 
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TABLE 1 

Log of interferometric observations for 70 Vir 


UT Date 

Baseline 

# of Obs 

Calibrators 

2013/04/03 

Sl/El 

5 

HD 119288, HD 121560, HD 122386 

2013/04/04 

Sl/El 

4 

HD 119550, HD 113022 

2014/04/20 

Wl/El 

4 

HD 119288, HD 119550 


Note. — Calibrator angular diameters from JSDCQ are: 0HD113O22 = 
0.39±0.03, flHDii9288 = 0.39±0.03, (9hdii9550 = 0.38±0.03, 6hdi21560 = 
0.46 ± 0.03, and 0HD122386 = 0.49 ± 0.03 mas. For details on the interfer¬ 
ometric observations, see m 

'‘Available at http://www.jmmc.fr/searchcal 


system for both the star and the known planet. Sec¬ 
tion [2] describes new interferometric observations ob¬ 
tained using the CHARA Array. Section [3] combines 
these measurements with a new spectroscopic analysis 
of Keck/HIRES data to determine fundamental stellar 
properties of 70 Vir. Section S] presents the addition 
of ^ 60 Keck/HIRES RV measurements to the existing 
time series, a revised Keplerian orbital solution, and the 
calculation of an accurate transit ephemeris. Section [5] 
uses the greatly improved system parameters to calcu¬ 
late the extent of the HZ and discusses the prospect of 
HZ planets in the system. Section [5] describes 19 years of 
precision robotic photometry that both rule out a plane¬ 
tary transit and show that the long-term stellar activity 
is constant within 0.004 mag. We provide concluding 
remarks in Section [3 

2. INTERFEROMETRIC OBSERVATIONS 

70 Vir (HD 117176; HR 5072; HIP 65721 ) is a bright 
{y = 4.97; H = 3.24: [Johnson et ^11968) and nearby 
(Hivv arcos parallax of 55.60 ± 0.24 mas: Ivan Leeiiwi^ 
120071 1 star. Our interferometric observations of 70 Vir 
were conducted at the Georgia State University’s Center 
for High Angular Resolution Astronomy (CHARA) Ar¬ 
ray and the Classic beam combiner in two-telescope mode 
operating in iJ-band (centra l wavelength Ac = 1.67^m) 
(jten Briimmelaar et al.l[2?i?)^ . We collected a total of 13 
observations over the course of three nights: two nights 
in April 2013 using the SlEl pair of telescopes, and one 
night in April 2014 using the WlEl pair of telescopes. 
The SlEl and ElWl are the two longest telescope con¬ 
figurations available at CHARA, with baselines B (dis¬ 
tances between two telescopes) of Hsiei = 330 m and 
JJwiEi = 313 m. 

Calibrator stars were observed in bracketed se- 
quences with 70 Vir We use the SearchCal software 
(|Bonneau et al.l 120061 1201111 to select calibrator stars 
based on their proximity in the sky with respect to 
the science target (within ~ 8 degrees). We observe 
a total of hve calibrators with estimated angular sizes 
Best < 0.5 mas in order to minimize systematic errors 
whic h could be introduced by the calibrator’s estimated 
sizes ()van Belle fc van Bellell200^ . A log of the observa¬ 
tions along with calibrator information can be found in 
Table m 

Calibrated data are used to determine the stellar uni¬ 
form disk angular diameter 0 ud and limb darkened an- 
gu lar diameter 6*T,n by h t ting the functions expressed 
in iHanburv Brown et al.l (1197411 . An estimate of the 
star’s temperature and gravity based on spectra are used 
to determine 77-band limb-darkening coefficients from 
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Fig. 1.— Plot of calibrated interferometric visibilities and the 
limb-darkened angular diameter fit for 70 Vir. The interferometric 
observations are described in Section [2] 

iClaret fc Bloem^ (|2011ll . These quantities are iterated 
upon with the final stellar parameters (see Section [3]) to 
determine the final coefficient used in th e limb darkened 
dianieter so lution, = 0.3512 (e.g., see lBovaiian et all 
We measure the angular diameter of 70 Vir 
to be 6 »ud = 0.967 ± 0.004 and 6 »ld = 0.998 ± 0.005 mil- 
liarcseconds. Figure [I] shows the interferometric data 
along with the best fit limb-darkened visibility function; 
the angular diameter measurements may be found in Ta¬ 
ble 1 

3. STELLAR PROPERTIES 
3.1. Interferometry 

The stellar angular diameter measured with interfer¬ 
ometry in Section [3] may be used in combination with 
the trigonometric parallax from Hippareos to derive the 
physical linear radius of the star, R, using trigonome¬ 
try. The most direct way to measure the effective tem¬ 
perature of a star is via the Stefan-Boltzmann equation, 
L = 47r7?^ crUgg, rearranged to yield 

Teff = 2341(^’bol/0L)°■"^ (1) 

where the constant 2341 absorbs the conversion constants 
assuming using units of Fboi in 10“® erg s“^ cm“^ and 
the limb-darkened angular diameter 0 ld in milliarcsec- 
onds (mas). 

The bolometric flux is measu red by norma lizing a 
G5 V spectral template from the iPicklesI l|1998ll library 
to broad-band photometry and spectrophotometry in 
the literature. Details of this method are described in 
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TABLE 2 

Stellar properties for 70 Vir 


Parameter 

Value 

Spectroscopic 

Value 

Interferometric Section Reference 

©uD (mas) . 


0.967 ±0.004 





Sld (mas) . 


0.998 ± 0.005 





G)oI (10“® erg cm” 

-2) 

28.050 ±0.562 

m \ 

Luminosity (Lq) . 


2.827 ±0.062 

m \ 

Radius R* (R©) . 

1.94 ±0.05 

1.9425 ± 0.0272 

m 3 \ I 

FT] 

Teff (K) . 

5439 ± 44 

5393 ± 30 

[ 



[Fe/H] . 

.. -0.09 ±0.03 



3.2 


sin i (km s . 

1.56 ±0.50 


[ 

3.2 


logs . 

3.90 ±0.06 



3.2 


Mass (M0) . 

1.09 ±0.02 



3.2 


Age (Gyr) . 

7.77 ±0.51 



3.2 



Note. — For details, see i[3l 


Ivan Belle et al.l (|2008ll : Ivon Braun et al.l (1201411 . For 


Johnson & MorganI ( 

195311: IGutierrez-Moreno et alJ 

(196611; Johnson et al. 

(1966H: Gowlev et al.l (1196711: 

Pfleiderer et al.l (I1966H: iJerzvkiewicz & Serkowski 

(196611: 

Piirola (1197611: 1 Johnson & Harris (1195411: 

Areue 

(HH; 

Serkowskil (119611: HiiggkvLst & Oial (11966 

1: Oia 

(198,511; IMermilliodI (I1986H: iJennens & 1 

lelfei ( 

1975): 

Moffett & BarnesI (1197911: iDucatil (1200211: 

Johnson et alJ 

(1196811: iBeichman et al.l (1198811: Gutri et al.l (200311: 

Gezari et al.l (1199911: Oial (119961): iDeanl (1198111: Olsen 

(11994): 

(TT^): 

Jasevicius et al.l (1199011: Rufener & Nicolet 

HiiggkvLst & Oial (1197011: iKornilov et al.l (199111: 


use th e spectrophotome t ry da t a from iKh aritono y et al.l 
(|1988ll : iGlushneva et al.l (1199811 1 iBurnashevI ( 198511 . 


The fit (see Figure[2]) produces a bolometric flux Fboi = 
(28.050 ± 0.0248) x 10“® erg s“^ cm“^. We note that 
the quoted uncertainty is statistical only and thus does 
not account for absolute errors in the templates, uncer¬ 
tain photometric zero-points, o r other effects such as th e 
ones outlined in Section 2.2 of Ivon Braun et HI (j2014[ l. 
We follow th e argu ments in Sections 3.2.1 and 3.2.2 of 
iBohlin et al.l (1201^ and add a 2% error in quadrature 
to account for a more realistic represe ntation of the true 
uncer tainties (see also the appendix in iBessell fc Murph^ 
l2012tl . The final Fboi and associated uncertainty values 
are presented in Table [21 along with the derived from 
Equation [T] using Fboi and 0ld- 


3.2. Spectroscopy 

We model two spectra of 70 Vir, taken on 3 
July 2009 with Keck/HIRES. The spectra are mod¬ 
eled using the Spectroscopy Made Easy (SME) pack - 
age (IValenti fc Piskuno^ll990l: IValenti fc Fischer! 1200511 . 
SME employs an iterative mode using results of the 
model atmosphere analysis i n combination with the 
Yonsei-Yale model isochrones (|Demaraue et al.l 1200411 in 
order to produce self- consistent results w ith the mea¬ 
sured surface gravity (IValenti et al.ll 200 '^ . The results 
of the spectroscopic modeling (surface gravity log 5 , ro¬ 
tational velocity ?;sini, atmospheric abundance [Fe/H], 
and effective temperature Teff) and stellar isochrone so¬ 
lution (mass M*, radius i?*, and age) are presented in 
Table |21 The values for effective temperature measured 
with interferometry and derived with spectroscopic mod¬ 
eling agree well. Likewise, there is excellent agreement 



Wavelength (/im) 

Fig. 2.— Upper panel. ( 1) th e blue curve is a G5 V spectral 
template from the IPicklesI II1998I 1 library, (2) the red crosses are 
literature photometry and spectrophotometry data of 70 Vir with 
errors in y-direction and filter bandwidths in x-direction, (3) the 
black X-shapes represent the specific flux values of the spectral 
template at the central wavelength of the filter of the respective 
literature photometry data point. Lower panel: the red crosses 
represent the fractional residuals to the fit. The high density of 
photometry points toward the blue end of the spectrum stems from 
the three spectrophotometry data sets in the literature for this star. 
For more details, see Section l3.ll 

with the radius predicted by model isochrones and the 
directly measured radius with interferometry. The intro¬ 
duction of the interferometric data set not only allows for 
an empirically based consistency check with the results 
from stellar atmosphere and evolutionary codes, but also 
reduces the uncertainties in the stellar parameters be¬ 
yond the capabilities of current methods employing mod¬ 
els (crTeff and cri?* are 85% and 47% lower respectively). 

3.3. Stellar Abundances 

There are at least 17 different groups who have mea- 
sured the s tellar abu ndances in 70 Vir, for example 
I Zhao et ^ (1200211 and Ida Silva et all (1201111 . Given the 
close proximity of the host star (~i 8 pc) and bright 
V magnitude, elements from lithium to europium have 
been measu r ed wi thin 70 Vir. Per the analysis of 
iHinkel et al.l (I2014I1 . the abundance measurements as 
determined by each group were renormalized to the 
iLodders et al.l (|2009ll solar abundance scale. The max¬ 
imum variation between datasets for each element, or 
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the spread, was determined in order to better charac¬ 
terize the consistent measurement of that element abun¬ 
dance. For 70 Vir, the median value for [Fe/H] = -0.01 
dex. while the sp read was 0.41 dex, since the renormal¬ 
ized [Lair^ (1198511 determination foun d [Fe/H ] = -0.2 dex 
and renormalized I Valenti &: Fischeil (1200811 determined 
[Fe/ H] = 0.2 dex. In t he analysis of the Hypatia Cata¬ 
log (|Hinkel et al.ll20l3l . to ensure abundances were co- 
pacetic, a star was not considered when the spread of 
the catalog abundances was greater than the error bar 
associated with that element, in this case ± 0.05 dex. 
Therefore, because of the discrepancy in [Fe/H] between 
datasets, 70 Vir was not incorporated into the analysis 
of Hypatia. 

Despite the ~30 elements that have been measured 
in 70 Vir, the majority of them suffer from inconsistent 
measurements between groups. The only elements that 
do not have a spread greater than the respective error 
bar are for cases where only one literature source has 
measured that element, or when the spread is 0.0 dex. 
In these cases, the renormalized abun dances are [Li/Fe] 
= -1 54 dex and [Eu/Fe] = 0.01 dex (IGonzalez fc La\^ 
UMl); [ K/Fe] = -0.21 dex (IZbao et 31.11200211 : [Zrll/Fe] 
= -0.13 (iMashonkina fc GehrenI 1200011: [VIl /Fel = -0.09 
dex and [Crll/Fe] = -0.42 dex (lTake'^1200711 : and [Sr/Fe] 
= 0.04 dex, fZr/Fe] = 0.09 dex, [Cell/Fe] = 0.08 dex 
(|da Silva et al.llMTlIl . While many of these abundances 
are sub-solar, not much can be said given their varying 
nucleosynthetic origins. 

4. A REFINED PLANETARY ORBIT 

Here we present new RV data for 70 Vir, a revised 
Keplerian orbital solution for the planet, and an accurate 
transit ephemeris for 70 Vir. 

4.1. Spectra Acquisition 

Previously published data for 70 Vir includes 169 mea¬ 
surements acquired with the Hamilton Echelle Spec¬ 
trograph (|Vogtl[l98^ on the 3.0m Shane Telescope at 
Lick Observatory (iMarcv &: Butleil Il996t iButler et al.l 
120061: lEischer et al.ll2ni4l l and 35 measurements acqu ired 
with the ELODIE spectrograph (jBaranne et al.l[T996ll on 
the 1.93m telesco pe at Observatoire de Haute-Provence 
(|Naef et al.l 1200411 . We add to these time series 59 new 
measurem ents acquired wi th the HIRES echelle spec¬ 
trometer (|Vogt et al.11199411 on the 10.0m Keck 1 tele¬ 
scope. H&K emission measured from the Keck spectra 
show that 70 Vir is a relatively quiet star. We show the 
complete dataset of 263 measurements in Table 13 where 
the fourth column indicates the source of the measure¬ 
ments. These data represent a baseline of ^ 26 years of 
monitoring 70 Vir. 


TABLE 3 

70 Vir Radial Velocities 


Date 

(JD - 2,440,000) 

RV 

(7 

(ms“^) 

Telescope 

Instrument 

7195.02270 

-155 

13 

Hamilton 

7195.02790 

-146 

13 

Hamilton 

7195.03310 

-146 

13 

Hamilton 

7224.84280 

145 

10 

Hamilton 

7224.84810 

146.0 

9.6 

Hamilton 

7224.85300 

144 

10 

Hamilton 

7373.71740 

109.7 

9.4 

Hamilton 


TABLE 3 — Continued 


Date 

(JD - 2,440,000) 

RV 

(ms-l) 

a 

(ms“^) 

Telescope 

Instrument 

7373.72130 

115.1 

9.5 

Hamilton 

7373.72640 

114.7 

8.5 

Hamilton 

7578.89680 

232.3 

9.3 

Hamilton 

7710.75280 

429 

11 

Hamilton 

7965.00190 

4 

12 

Hamilton 

8018.84930 

-130 

13 

Hamilton 

8376.80300 

-66 

12 

Hamilton 

8437.74980 

-55 

12 

Hamilton 

8670.98610 

-73 

11 

Hamilton 

8670.99340 

-82 

14 

Hamilton 

8746.68150 

254 

14 

Hamilton 

8746.70340 

291 

13 

Hamilton 

8781.77040 

-12 

11 

Hamilton 

8781.79240 

4 

10 

Hamilton 

8847.67430 

-37 

11 

Hamilton 

8993.08360 

442.3 

9.9 

Hamilton 

9068.91280 

-110 

10 

Hamilton 

9068.93410 

-110.3 

9.8 

Hamilton 

9096.85900 

263 

10 

Hamilton 

9096.88040 

274 

11 

Hamilton 

9124.77090 

94 

10 

Hamilton 

9124.79250 

98 

11 

Hamilton 

9172.70420 

-141 

11 

Hamilton 

9349.09110 

305 

11 

Hamilton 

9411.97090 

-128 

11 

Hamilton 

9469.81290 

216 

13 

Hamilton 

9767.00320 

-120.8 

5.2 

Hamilton 

9768.97949 

-119.3 

5.6 

Hamilton 

9793.94141 

176.7 

4.7 

Hamilton 

9793.95117 

187.2 

4.2 

Hamilton 

9793.96191 

185.4 

4.2 

Hamilton 

9793.97168 

188.2 

3.9 

Hamilton 

10087.93164 

-162.9 

4.4 

Hamilton 

10089.01953 

-162.8 

4.5 

Hamilton 

10089.97363 

-157.4 

3.7 

Hamilton 

10091.08594 

-163.3 

4.4 

Hamilton 

10121.02539 

-113.9 

5.5 

Hamilton 

10121.04785 

-105.5 

6.1 

Hamilton 

10122.01562 

-108.8 

7.3 

Hamilton 

10122.03711 

-115.9 

6.3 

Hamilton 

10124.98535 

-83.0 

4.6 

Hamilton 

10125.00977 

-92.9 

4.4 

Hamilton 

10125.86133 

-84.6 

5.4 

Hamilton 

10127.01758 

-77.8 

6.5 

Hamilton 

10127.06055 

-78.8 

6.9 

Hamilton 

10129.03223 

-52.2 

4.6 

Hamilton 

10129.05469 

-56.4 

4.8 

Hamilton 

10144.95215 

206.6 

3.8 

Hamilton 

10145.94727 

231.7 

5.8 

Hamilton 

10145.95996 

242 

10 

Hamilton 

10148.95410 

320.3 

5.6 

Hamilton 

10148.98731 

316.2 

5.8 

Hamilton 

10150.92188 

386.5 

4.8 

Hamilton 

10150.94434 

367.0 

4.7 

Hamilton 

10172.90918 

133.3 

4.6 

Hamilton 

10172.92969 

139.5 

8.7 

Hamilton 

10173.88769 

111.6 

5.1 

Hamilton 

10179.78027 

23.8 

4.1 

Hamilton 

10180.76367 

7.5 

3.8 

Hamilton 

10180.76953 

4.8 

4.3 

Hamilton 

10186.82031 

-41.7 

5.3 

Hamilton 

10186.84375 

-52.4 

5.9 

Hamilton 

10199.80957 

-132.2 

5.0 

Hamilton 

10199.83203 

-132.4 

4.9 

Hamilton 

10200.82910 

-142.9 

4.8 

Hamilton 

10200.85059 

-143.7 

5.4 

Hamilton 

10201.81641 

-135.7 

5.5 

Hamilton 

10201.83887 

-142.2 

5.9 

Hamilton 

10202.81836 

-144.4 

5.3 

Hamilton 

10202.84082 

-152.1 

4.7 

Hamilton 

10214.76856 

-172.6 

5.0 

Hamilton 

10214.79102 

-171.4 

5.1 

Hamilton 

10234.79102 

-128.4 

5.7 

Hamilton 

10234.81250 

-118.0 

5.5 

Hamilton 

10263.66894 

261.4 

3.4 

Hamilton 

10502.97852 

412.0 

4.7 

Hamilton 
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TABLE 3 — Continued 


TABLE 3 — Continued 


Date 

(JD - 2,440,000) 

RV 

(ms“^) 

(T 

(m s~^) 

Telescope 

Instrument 

Date 

(JD - 2,440,000) 

RV 

(ms“^) 

(T 

(ms~^) 

Telescope 

Instrument 

10537.88281 

-72.0 

4.4 

Hamilton 

14253.74064 

198.9 

5.0 

Hamilton 

10563.75195 

-154.4 

4.8 

Hamilton 

14254.73845 

183.6 

4.9 

Hamilton 

10614.75391 

288.4 

5.2 

Hamilton 

14547.96363 

-155.9 

4.9 

Hamilton 

10614.75879 

283.9 

4.7 

Hamilton 

14574.87404 

98.5 

4.8 

Hamilton 

10655.68848 

-80.3 

4.9 

Hamilton 

14864.06816 

-135.3 

6.9 

Hamilton 

10794.06055 

-162.7 

3.8 

Hamilton 

14864.07598 

-131.4 

7.6 

Hamilton 

10978.75684 

395.0 

5.8 

Hamilton 

15229.03541 

-165.0 

7.2 

Hamilton 

10978.77148 

408.4 

5.2 

Hamilton 

15229.03803 

-165.1 

7.3 

Hamilton 

11005.70801 

-81.5 

5.8 

Hamilton 

15312.79659 

46.8 

4.4 

Hamilton 

11207.01269 

466.0 

5.2 

Hamilton 

15312.79906 

47.0 

4.2 

Hamilton 

11303.82129 

43.7 

5.6 

Hamilton 

10150.56740 

399.0 

7.0 

ELODIE 

11700.72559 

-41.3 

9.5 

Hamilton 

10207.47330 

-124.0 

7.0 

ELODIE 

11700.74805 

-31.2 

6.8 

Hamilton 

10212.49560 

-121.0 

7.0 

ELODIE 

11946.02832 

-124.0 

4.7 

Hamilton 

10263.35220 

295.0 

7.0 

ELODIE 

11969.94727 

-157.3 

5.8 

Hamilton 

10267.36290 

411.0 

8.0 

ELODIE 

12041.82715 

94.9 

6.0 

Hamilton 

10477.67260 

-19.0 

7.0 

ELODIE 

12041.83301 

90.6 

5.4 

Hamilton 

10532.52620 

28.0 

7.0 

ELODIE 

12055.77441 

-82.7 

5.4 

Hamilton 

10535.54120 

-3.0 

7.0 

ELODIE 

12055.80371 

-76.9 

5.0 

Hamilton 

10557.54150 

-104.0 

7.0 

ELODIE 

12072.70117 

-164.1 

5.9 

Hamilton 

10561.51440 

-129.0 

7.0 

ELODIE 

12072.73047 

-153.7 

6.0 

Hamilton 

10582.41520 

-90.0 

7.0 

ELODIE 

12122.69629 

71.7 

4.8 

Hamilton 

10587.41230 

-62.0 

7.0 

ELODIE 

12333.97461 

-122.3 

6.0 

Hamilton 

10623.40340 

506.0 

7.0 

ELODIE 

12334.95410 

-121.8 

6.9 

Hamilton 

10858.59380 

515.0 

7.0 

ELODIE 

12335.96387 

-123.3 

4.8 

Hamilton 

10885.57180 

0.0 

8.0 

ELODIE 

12338.95898 

-99.6 

5.7 

Hamilton 

10940.42360 

-52.0 

7.0 

ELODIE 

12345.93164 

-55.4 

5.7 

Hamilton 

10967.38640 

395.0 

7.0 

ELODIE 

12345.95606 

-53.0 

6.3 

Hamilton 

11024.34770 

-121.0 

8.0 

ELODIE 

12348.95215 

-25.0 

6.4 

Hamilton 

11239.62410 

-32.0 

7.0 

ELODIE 

12348.97656 

-29.9 

4.9 

Hamilton 

11300.41400 

38.0 

7.0 

ELODIE 

12427.76367 

-176.5 

5.2 

Hamilton 

11325.41920 

505.0 

7.0 

ELODIE 

12428.80859 

-178.1 

4.3 

Hamilton 

11588.57810 

-20.0 

7.0 

ELODIE 

12449.72754 

-141.5 

6.1 

Hamilton 

11623.51300 

-103.0 

7.0 

ELODIE 

12449.73047 

-134.8 

5.8 

Hamilton 

11653.51820 

88.0 

7.0 

ELODIE 

12723.92578 

465.2 

5.6 

Hamilton 

11692.42330 

135.0 

7.0 

ELODIE 

12723.94922 

459.0 

6.1 

Hamilton 

11956.64950 

-94.0 

7.0 

ELODIE 

12796.69531 

-154.1 

4.3 

Hamilton 

12038.51990 

208.0 

7.0 

ELODIE 

12796.71777 

-152.5 

4.5 

Hamilton 

12356.57740 

133.0 

7.0 

ELODIE 

12796.74023 

-147.3 

4.0 

Hamilton 

12361.58960 

246.0 

7.0 

ELODIE 

13021.07031 

-166.8 

5.2 

Hamilton 

12413.50440 

-73.0 

8.0 

ELODIE 

13022.02494 

-166.9 

5.0 

Hamilton 

12719.55560 

446.0 

7.0 

ELODIE 

13069.87695 

400.5 

4.3 

Hamilton 

12748.52720 

50.0 

7.0 

ELODIE 

13080.93945 

367.0 

4.7 

Hamilton 

12751.53870 

10.0 

7.0 

ELODIE 

13080.96191 

360.2 

4.2 

Hamilton 

12772.41720 

-103.0 

7.0 

ELODIE 

13102.86621 

-51.7 

4.8 

Hamilton 

12776.49050 

-105.0 

7.0 

ELODIE 

13130.82324 

-171.7 

4.3 

Hamilton 

13933.82196 

-217.7 

1.0 

HIRES 

13130.82617 

-166.9 

5.1 

Hamilton 

13933.82264 

-221.0 

1.0 

HIRES 

13162.69434 

-72.8 

4.5 

Hamilton 

13933.82336 

-217.8 

1.0 

HIRES 

13162.71875 

-61.4 

4.9 

Hamilton 

14085.16632 

-203.3 

1.3 

HIRES 

13363.05150 

-165.9 

4.7 

Hamilton 

14085.16702 

-204.2 

1.3 

HIRES 

13389.06953 

-111.7 

4.7 

Hamilton 

14085.16775 

-203.6 

1.2 

HIRES 

13391.00264 

-106.2 

5.0 

Hamilton 

14130.12242 

321.0 

1.1 

HIRES 

13392.04633 

-97.8 

4.6 

Hamilton 

14130.12308 

318.4 

1.1 

HIRES 

13393.02009 

-92.3 

4.7 

Hamilton 

14130.12372 

319.2 

1.1 

HIRES 

13403.04659 

8.6 

4.6 

Hamilton 

14131.11653 

290.6 

1.2 

HIRES 

13436.90184 

205.2 

5.2 

Hamilton 

14131.11724 

291.0 

1.2 

HIRES 

13438.93611 

150.6 

5.1 

Hamilton 

14131.11797 

291.7 

1.1 

HIRES 

13440.88377 

122.2 

5.5 

Hamilton 

14139.10185 

90.7 

1.2 

HIRES 

13441.98706 

100.5 

4.8 

Hamilton 

14139.10250 

91.1 

1.2 

HIRES 

13475.82543 

-145.5 

4.8 

Hamilton 

14139.10312 

87.7 

1.2 

HIRES 

13476.79240 

-152.8 

5.1 

Hamilton 

14641.84867 

-237.2 

1.2 

HIRES 

13477.81099 

-161.7 

4.9 

Hamilton 

14807.16184 

12.3 

1.6 

HIRES 

13478.84384 

-162.9 

4.9 

Hamilton 

14809.16606 

46.7 

1.3 

HIRES 

13479.85233 

-159.4 

4.9 

Hamilton 

14811.16959 

95.0 

1.2 

HIRES 

13501.82814 

-138.0 

5.1 

Hamilton 

14865.09274 

-199.3 

1.2 

HIRES 

13565.71260 

-17.6 

5.5 

Hamilton 

14866.06558 

-209.0 

1.3 

HIRES 

13753.02186 

24.3 

7.3 

Hamilton 

14868.10796 

-213.6 

1.5 

HIRES 

13756.07980 

66.1 

5.2 

Hamilton 

14927.00429 

68.3 

1.4 

HIRES 

13773.99421 

448.6 

4.9 

Hamilton 

14963.76892 

-54.2 

1.4 

HIRES 

13843.87891 

-148.6 

4.4 

Hamilton 

14983.92452 

-212.9 

1.4 

HIRES 

14134.03032 

276.3 

4.5 

Hamilton 

14984.93831 

-211.1 

1.2 

HIRES 

14169.93907 

-150.6 

4.8 

Hamilton 

14985.93592 

-217.3 

1.2 

HIRES 

14196.85956 

-149.1 

4.6 

Hamilton 

14986.92704 

-216.2 

1.2 

HIRES 

14219.84280 

24.7 

4.9 

Hamilton 

14987.92941 

-221.9 

1.2 

HIRES 

14253.73413 

199.8 

5.2 

Hamilton 

14988.87370 

-227.3 

1.4 

HIRES 

14253.73740 

191.8 

4.9 

Hamilton 

15016.74253 

-202.7 

1.3 

HIRES 
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TABLE 3 — Continued 


Date 

(JD - 2,440,000) 

RV 

(ms-l) 

(J 

(m s“^) 

Telescope 

Instrument 

15042.73366 

58.4 

1.2 

HIRES 

15043.74347 

74.7 

1.3 

HIRES 

15044.77665 

100.2 

1.2 

HIRES 

15256.94773 

-182.2 

1.3 

HIRES 

15286.01509 

301.1 

1.4 

HIRES 

15311.80566 

-30.5 

1.3 

HIRES 

15312.79832 

-39.2 

1.1 

HIRES 

15314.85494 

-74.7 

1.2 

HIRES 

15375.75214 

-167.0 

1.2 

HIRES 

15405.73625 

377.3 

1.1 

HIRES 

15585.18842 

-241.5 

1.2 

HIRES 

15607.18024 

-180.8 

1.3 

HIRES 

15636.97057 

333.8 

1.4 

HIRES 

15636.97113 

330.3 

1.3 

HIRES 

15636.97172 

332.3 

1.2 

HIRES 

15636.97232 

329.7 

1.2 

HIRES 

15706.73340 

-238.4 

1.0 

HIRES 

15707.73062 

-237.0 

1.2 

HIRES 

15731.94048 

-113.1 

1.3 

HIRES 

15792.72022 

-177.0 

1.3 

HIRES 

15792.72080 

-170.0 

1.3 

HIRES 

15792.72139 

-172.8 

1.2 

HIRES 

15993.02461 

379.9 

1.3 

HIRES 

15993.02698 

382.7 

1.4 

HIRES 

15993.03026 

386.1 

1.5 

HIRES 

16145.72624 

-185.9 

1.2 

HIRES 

16488.73474 

-142.7 

1.1 

HIRES 

16675.18710 

24.3 

1.4 

HIRES 


4.2. Keplerian Orbital Solution 

The revised Keplerian orbital solution to the RV data 
in Table |3] used RVLIN; a partiall y linearized, least- 
square s fitting procedure described in iWright &: HowardI 
(|2009ll . Parameter uncertainties were estimated using 
the BOOTTRAN bootstrapping routines described in 
iWang et al.l (|2012ll . The resulting orbital solution is 
shown in Table 2] and in Figure [3] 


TABLE 4 

Keplerian Orbital Model 


Parameter 

Value 

70 Vir b 

P (days) 

116.6926 ±0.0014 

Tc “ (JD - 2,440,000) 

16940.258 ± 0.084 

Tp >> (JD - 2,440,000) 

7239.7091 ± 0.11 

e 

0.399 ± 0.002 

OJ (deg) 

358.8 ±0.3 

K (ms-l) 

315.7 ±0.7 

Mp sin 2 {Mj) 

7.40 ±0.02 

a (AU) 

0.481 ±0.003 

System Properties 

7 (ms 

22.94 ±0.59 

Measurements and Model 

t^ohs 

263 

rms {ms“^) 

6.08 

X?ed 

1.16 


^ Time of mid-transit. 

Time of periastron passage. 

For each bootstrapping realization, the fit produces off¬ 
sets for each dataset with respect to the Lick Hamilton 
dataset. These are fit as two additional free parameters 
in the Keplerian orbital fit described above. We find the 


offsets to be 48.4 and -74.6 ms“^ for data from ELODIE 
and HIRES respectively. The and rms scatter of the 
residuals (see Table a are consistent with the measure¬ 
ment uncertainties shown in Tabled Note that we added 
a stellar jitter noise component of 3 ms~^ in quadratur e 
with the measurement uncertainties (jBiitler et al.ll200^ . 
We hnd no evidence for a linear RV trend in the fit resid¬ 
uals shown in the right panel of Figure [3l 

We performed a further analysis of the data to search 
for signatures of possible additional planets. Figure a 
shows the Lomb-Sc argle periodogram (|Horne fc Baliunaa 
119861 : [^arglelll98^ of the best-fit residuals, which shows 
no dominant peak. Since the lower-precision ELODIE 
and Hamilton data might obscure a low-amplitude signal 
detectable in the HIRES data, we have also examined the 
HIRES data alone. Figure |4] shows the periodogram of 
the residuals of the HIRES data to the best ht shown in 
Table 01 

The presence of many peaks of similar amplitude in 
these periodograms is consistent with there being many, 
low-mass planets of similar RV semi-amplitude in the 
data, but also consistent with noise. Since the rms of 
the residuals is consistent with both expectations and 
measurements of the uncertainties, there is no reason to 
expect the former, so we conclude that these data contain 
no evidence of additional periodic astrophysical signals. 

4.3. Transit Ephemeris Refinement 

The transit mid-point epoch shown in Table 0] was cal¬ 
culated with a Monte-Carlo bootstrap, which propagates 
the uncertainty in this orbital parameter to the time of 
the transit. This method produces the most accurate 
ephemeris since the transit times are calculated as part 
of the orbital fit. Note that if the planet does not tran¬ 
sit then the transit mid-point epoch may be considered 
the time of inferior conjunction. The predicted transit 
properties of a system depend sensitively on the stellar 
radius as well as the planetary parameters. We adopt 
our interferometric measured radius from Table [2] which 
has an uncertainty of only 1.5%. The minimum mass of 
the planet is larger than a Jupiter mass and we approxi¬ 
mate the radius of the planet as Rp = l.O R.j, based upon 
the ma ss-radius relationship described bv iKane fc Gelinol 
(|2012all . These properties, combined with the orbital so¬ 
lution of Table 01 result in a transit probability of 2.27%, 
a predicted transit duration of 0.66 days, and a tran¬ 
sit depth of 0.3%. The transit mid-point uncertainty 
shown in Table0]is 0.084 days, or 121 minutes. Therefore 
the transit window is dominated by the transit duration 
rather than the mid-point uncertainty, which is a favor¬ 
able scenario for photometric follow-up. Our procedure 
is to use a calculated value for Tc as close as possible 
to the conclusion of observations. However, the baseline 
of the RV observations described in Section 14.11 is long 
enough such that there is very little increase in the size 
of this transit window for the foreseeable future. The 
uncertainty in the predicted transit time subsequent to 
that shown in Table 0] has an uncertainty that is less 
than a minute larger. iKane fc von Brauiil (|2008f l have 
also shown that the transit probability is a strong func¬ 
tion of both the eccentricity and the argument of pe¬ 
riastron. For example, if the eccentricity of the planet 
were zero, the transit probability and duration would be 
1.92% and 0.71 days respectively. Thus the orientation 
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Fig. 3.— Left All 263 RV measurements from three different instruments (see Table 0 for 70 Vir phased on the new orbital solution 
shown in Table [4| RV offsets between datasets have been accounted for in this figure. Right Residual velocities with respect to the best 
orbital solution. 
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Fig. 4. — Left: A Lomb-Scargle periodogram of the residuals to the Keplerian orbital fit shown in Table [JJ using all of the available data 
from Table 0 The fourier powers shown in the figure are consistent with noise in the data. Right: The same analysis repeated using only 
the HIRES data, showing that our data do not reveal the presence of additional planets. 


of the 70 Vir b orbit results in a slightly enhanced transit 
probability relative to a circular orbit. 

5. SYSTEM HABITABLE ZONE 

The fundamental stellar parameters from Table [2] pro¬ 
vide the means to investigate the HZ of the system and 
the potential for terrestrial planets in that re gion. Previ¬ 
ous st udies of the 70 Vir HZ include those of ijones et al.l 
(|2006fl . who calculated HZ boundaries for a selection of 
known exoplanetary systems using the esti mated ages of 
stars to determine on-going habitability. iSandor et al.l 
()2007fl studied stability regions in the 70 Vir system and 
concluded that the system is unlikely to host HZ plan¬ 
ets. Th ese previous s t udies used the older HZ bound¬ 
aries of iKasting et al.l (|1993f) . Here we revisit the HZ 
properties of 70 Vir using the revised system parameters 
pr esented here along with the u pdated HZ calculations 
of lKoDparapu et ^ (1201,1112014I 1. 

We adopt the definitions of “con servative” a n d “op - 
timistic” HZ models described by iKane et all (120131) . 
The conservative HZ use boundaries based upon run¬ 
away and maximum greenhouse climate models, whereas 
the optimistic HZ extends these boundaries based on as¬ 


sumptions regarding the amount of time that Venus and 
Mars were able to retai n liquid water on their surfaces 
(iKopparapu et al.ll2013f l. The accuracy to which these 
boundaries can be determined re ly on robus t determi¬ 
nations of the stellar parameters (jKanell^Oldl) which, in 
this case, have exceptionally small related uncertainties 
(see Table [5]) such that the HZ boundary uncertainties 
are negligible. HZ calculations for all known exoplan¬ 
etary systems are available using the s ame methodol¬ 
ogy th rough the Habitable Zone Gallery (|Kane fc Gelinol 
1291 2aD . 

Figure [5] shows a top-down view of the 70 Vir system 
where the solid line indicates the Keplerian orbit of the 
planet using the orbital parameters of Table ID The HZ 
is depicted by the shaded region where the light gray 
represents the conservative HZ and the dark gray is the 
optimistic extension to the HZ. The conservative HZ cov¬ 
ers the region 1.63-2.92 AU from the host star and the 
optimistic HZ increases this region to 1.29-3.08 AU. 

Although the confirmed planet is clearly interior to 
the HZ, we performed stability simulations to investigate 
whether the relatively large mass of the planet and the 
eccentricity of its orbit exclude the presence of a hypo- 
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Fig. 5.— A top-down view of the 70 Vir system showing the ex¬ 
tent of the HZ calculated using the stellar parameters of Table [2] 
The conservative HZ is shown as light-gray and optimistic exten¬ 
sion to the HZ is shown as dark-gray. The revised Keplerian orbit 
of the planet from Table [4| is overlaid. 

thetical Earth-mass planet in the HZ. To accomplish this, 
we performed dynamical simulations using N- body inte¬ 
grati ons with the Mercury Integrator Package (jChambersI 
Il9^ . We adopted the hybrid symplectic/Bulirsch-Stoer 
integrator and used a Jacobi coordinate system, which 
provides more accurate results for multi -planet systems 


(|Wisdom h HolmanI 

199lHWisdoinll200^ except in cases 

of close encounters 

Chamber§||1999t'- We inserted the 


hypothetical planet in a circular orbit at each of the four 
optimistic and conservative HZ boundaries. The inte¬ 
grations were performed for a simulation of 10® years, in 
steps of 100 years, starting at the present epoch. 

Assuming that the system is coplanar with an incli¬ 
nation of 90°, our simulations show that the hypotheti¬ 
cal systems all remain stable for the full duration of the 
simulations. The eccentricity of the hypothetical planet 
oscillates over the course of the simulation with a range 
of 0.00-0.35, 0.04-0.30, 0.05-0.22, and 0.03-0.18 for the 
optimistic inner, conservative inner, conservative outer, 
and optimistic outer boundaries respectively. These ec¬ 
centricities do not necessarily rule out habitability of 
the planet, d epending on the dynamics of the planetary 
atmo sphere (iKane fc Gelind I2dl2bl : iWilliams fc Pollai^ 
l2002f) . Note that stability at the HZ boundaries does 
not guarantee stability within the boundaries as that is 
a complex function of orbital distance, phase, and eccen¬ 
tricity. 

Since the mass of the inner planet depends on the incli¬ 
nation of the system {MpSini = 7A0Mj), we performed 
simulations that determine the system inclination where 
the mass of the inner planet causes the orbit of the outer 
planet to become unstable. These stability threshold in¬ 
clinations for the four boundaries are 24°, 25°, 10°, and 
3° for the optimistic inner, conservative inner, conserva¬ 
tive outer, and optimistic outer boundaries respectively. 
Shown in Figure [6] are the simulation results at the sta¬ 


bility threshold inclination for the inner optimistic and 
conservative HZ boundaries. Each panel shows the ec¬ 
centricity oscillations for the 50,000 years leading up to 
the instability event. Even though the inner conservative 
HZ boundary is farther away from the inner planet, the 
orbital period at that boundary places the outer planet 
closer to an orbital resonance with the inner planet than 
an orbit at the inner optimistic HZ boundary. Thus the 
planet remains stable for less time at the former than the 
latter. 

6. PHOTOMETRIC OBSERVATIONS 

We have been monitoring 70 Vir for two decades 
with the T4 0.75 m automatic photoelectric telescope 
(APT) located at Fairborn Observatory in southern Ari¬ 
zona. The T4 APT observes in the Stromgren b and y 
pass bands with an EMI 9124QB photomultiplier tube 
(PMT) as the detector. The automated photometer has 
a Fabry lens placed behind the focal-plane diaphragm 
that projects a fixed image of the primary mirror (illumi¬ 
nated by the star) onto the photo-cathode of the PMT. 
Thus, slight motions of the star within the diaphragm 
during an integration do not translate into image motion 
on the PMT cathode. The instrumentation and observ¬ 
ing strategy result in the data being close to the pho¬ 
ton/scintillation noise limit with far less correlated noise 
than is typical of CCD photometry which suffers from, 
for example, intra-pixel sensitivity. The data are thus as¬ 
sumed to be uncorrelated in the subsequent analysis. As 
an additional verification of validity of this assumption, 
each of the APT integrations are divided into a series 
of 0.1 second integrations and saved for quality control 
and trouble shooting purposes. Histograms of the subin¬ 
terval data and computed Geneva statistics are used to 
verify that the data are Gaussian and determine if there 
are trends, cycles, spikes, or drops in the photon counts 
during the integration that require further investigation. 
The T4 APT, its photometer, observing techniques, data 
reduction procedures, and photometric p recision are de¬ 
scribed in further detail in iHenr^^ (jl999| l. 

The comparison star HD 117304 (Cl: V = 5.65, 
B — V = 1.05, KO HI) has been used for all 23 observ¬ 
ing seasons since 1993, while comparison star HD 112503 
(C2: V — 6.81, B — V = 0.47, F7 IV) has been used only 
for the past 19 observing seasons beginning in 1997 be¬ 
cause it was chosen to replace a previous comparison star 
recognized to be a low-amplitude variable after the first 
four years. T4 has acquired 2051 good differential ob¬ 
servations with Cl over the past 23 years and 1897 good 
observations with C2 in the past 19 observing seasons. 
During the course of our analysis, we recognized that 
comparison star Cl exhibited very low-amplitude vari¬ 
ability at times; therefore, in this paper, we present the 
results of our analysis of the 1897 differential magnitudes 
of 70 Vir with respect to comparison star HD 112503 
(C2). 

The 1897 differential magnitudes computed with C2 
are plotted in the top panel of Figure [71 To increase 
the precision of the differential magnitudes, we com¬ 
bined the b and y observations into a single {b + y)/2 
’’pass band.” We also normalized each observing season 
to have the same mean magnitude as the first, thus mak¬ 
ing our search for short-period variability and shallow 
transits more sensitive. The nightly normalized observa- 
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Inner Optimistic Habitable Zone Boundary 
Orbital Period = 512.5 days 
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Inner Conservative Habitable Zone Boundory 
Orbital Period = 727.9 days 
Orbital Inclination (Stability Threshold) = 25° 
Time Until Instability = 152800 years 
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Fig. 6.— Stability simulations for a hypothetical Earth-mass planet in the HZ of the 70 Vir system. Each panel shows the eccentricity 
oscillations for the 50,000 years leading up to the ejection of the outer planet. Top panel: The Earth-mass planet remains stable at the 
inner boundary of the optimistic HZ for system inclinations > 24°, otherwise the inner planet causes the outer planet to be ejected after 
350,000 years. Bottom panel: The Earth-mass planet remains stable at the inner boundary of the conservative HZ for system inclinations 
> 25°, otherwise the inner planet causes the outer planet to be ejected after ~ 150,000 years. 


tions scatter about their grand mean, indicated by the 
straight line in the top panel, with a standard devia¬ 
tion 0.00270 mag. This is slightly larger than our typi¬ 
cal measurement precision given above and may indicate 
slight residual variability in 70 Vir and/or the comp star 
HD 112503. 

The observations are replotted in the middle panel 
of Figure [71 where they have been phased with the 
time of conjunction and the orbital period from Table 
0) A least-squares sine fit on the 116.6926-day radial 
velocity period gives a formal semi-amplitude of just 
0.000352 ± 0.000076 mag, thus confirming that the pe¬ 
riodic radial velocity variations are due to planetary re¬ 
flex motion and not to intrinsic stellar brightness vari¬ 
ations (see. e.g..rQueloz et al.ll200ll : iPaulson et al.ll2?i?)3 : 
iBoisse et al.ll20lT ^ 

The observations within ±0.03 phase units of the pre¬ 
dicted transit time are plotted in the bottom panel of 
Figure 0 The solid curve shows the predicted transit 
phase (0.0), the transit depth (0.3% or 0.00325 mag), 
and transit duration (±0.003 phase units) computed as 
described in Section 14.31 above. The horizontal line be¬ 
low the transit window represents the transit mid-point 
uncertainty. While the second half of the transit window 
is not well covered by our observations, there are a total 
of 27 observations within the transit window that have a 


mean of —1.773035±0.000416 mag and 1870 observations 
that fall outside the transit window and have a mean of 
-1.772784±0.000063mag. Thus, our ’’observed” transit 
depth is —0.000251 ± 0.000421 mag, which is consistent 
with zero to three decimal places. Therefore, a central 
transit of the expected depth and duration occurring at 
the expected time can be ruled out at the 5cr level. Al¬ 
though the data sampling is sufficient to also constrain 
the absence of transits for almost all impact parame¬ 
ters, the number of data points within the corresponding 
transit durations will be less, thereby lessening the signif¬ 
icance of such constraints. For example, the largest gap 
in the photometry during the transit window corresponds 
to ~0.3 of the central transit duration. This means the 
impact parameter would need to be > 0.95 to have been 
completely missed by our data. Ruling out such a range 
of impact parameters would reduce the posterior transit 
probability from 2.27% to 0.11%. 

70 Vir is a magnetically inactive star with logR/j^ 
values of -4. 99 and -5.116 acco r ding t o iWright et al.l 
(1200411 and llsaacson fc Fischeil (120101 1. respectively. 
IWright et al.l ()2004| 1 give an estimated rotation period 
of 32 days for 70 Vir, based on the star’s activity level. 
However, no reliable rotation period for 70 Vir has 
been directly measured via rotational modulation of dark 
starspots or bright Ca H H and K regions across the face 
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Fig. 7.— Top: The 1997—2015 normalized differential observa¬ 
tions of 70 Vir with respect to comparison star C2, acquired with 
the T4 0.75 m APT between 1997 and 2015. Middle: The obser¬ 
vations phased with the 116.6926 day orbital period. The semi¬ 
amplitude of a least-squares sine fit to the phase observations is 
0.000140 di 0.000081 mag, consistent with the absense of light vari¬ 
ability on the radial velocity period and confirming planetary re¬ 
flex motion of the star as the cause of the radial velocity varia¬ 
tions. Bottom: The observations within ±0.03 phase units of the 
predicted transit time. The solid curve shows the predicted tran¬ 
sit time at phase 0.0, transit depth (0.3%) and duration (±0.003 
phase units) for a central transit of planet b. The short horizontal 
line segment represents the uncertainty in the time of tansit. Our 
photometry shows that central transits of the expected depth and 
duration likely do not occur. 

of th e star (e.g., iHenrv et al.lll997l l200Ct I Simpson et alJ 
l2010f) . We performed periodogram analyses of each indi¬ 
vidual observing season and of our data set as a whole, 
and, while we found suggestions of low-amplitude vari¬ 
ability, we could not identify any significant period that 
might be interpreted as the stellar rotation period. 

Finally, we look for long-term variability in 70 Vir. Un¬ 
fortunately, comparison star Cl has significant long-term 
variability of several mmag. The yearly-mean differential 
magnitudes (70 Vir - Cl) vary over a range 0.0076 mag 
and have a standard deviation of 0.0023 mag with respect 
to the grand mean. However, the yearly means of (70 Vir 
- C2) have both a smaller range and a smaller standard 
deviation, 0.0041 mag and 0.0012 mag, respectively (see 
Figure [8|). Without another good comparison star, we 
cannot determine whether the variability we see in Fig¬ 
ure [5] originates in 70 Vir, the comparison star, or a com¬ 
bination of both. Thus, we can only state that 70 Vir’s 
long-term variability has a range less than ^0.004 mag. 

7. CONCLUSIONS 

In an era where new planets are being regularly discov¬ 
ered via the transit method, the bright exoplanet host 
stars still largely belong to those planets which were dis¬ 



Observed Depth = -0.000251 +/— 0.000421 mag J 


C\2 -1.78 
% -1.77 
-1.76 


Fig. 8.— The 1997-2015 yearly mean differential magnitudes 
70 Vir minus comparison star C2. The error bars on the individual 
seasonal means are slightly smaller than the plotted points. The 19 
differential magnitudes vary over a total range of 0.0041 mag and 
have a standard deviation of 0.0012 mag with respect to the grand 
mean, indicated by the horizonal dotted line. Without another 
good comparison star, we are unable to determine whether the 
variability is intrinsic to 70 Vir or the comparison star, so these 
values can only be quoted as upper limits of long-term variability 
in 70 Vir. 

covered using the radial velocity method. These are sys¬ 
tems which thus provide the greatest access to follow-up 
investigations due to the relatively large signal-to-noise 
possibilities presented. Here we have presented new re¬ 
sults for the 70 Vir system which includes detailed char¬ 
acterization of the host star. Our direct measurements 
of the stellar radius show that, although slightly cooler, 
70 Vir is almost twice the size of the Sun. This is consis¬ 
tent with the star being older and more evolved than the 
Sun. Our new radial velocity data provide an improved 
Keplerian orbital solution for the planet and further evi¬ 
dence that there are unlikely to be further giants planets 
within the system. A terrestrial-mass planet may yet 
exist beneath our detection threshold and so, given the 
vastly improved stellar properties, we calculated the HZ 
boundaries and performed stability simulations within 
that region. Our simulations show that a terrestrial 
planet could remain in a stable orbit near the HZ in¬ 
ner edge for system inclinations >25° and close to the 
outer HZ edge for almost all system inclinations. Finally, 
our 19 years of APT photometry confirm that the star 
is quite stable over long time periods and there is no 
evidence that the b planet tran sits the host star (a “dis¬ 
positive null”, as described by IWang et al.l (|2012ll l. the 
timing of which we were able to accurately predict from 
our revised Keplerian orbit. The TERMS compilation 
of data for this system presented here means that it is 
now one of the better characterized systems in terms of 
stellar and planetary parameters. 
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